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Bounds  and  comparisons  for  various  performance  measures  in  tandem 
queueing  systems  and  loss  systems  are  studied. 


An  upper  bound  for  the  expected  delay  in  front  of  the  second  ser- 
ver is  found  for  a sequence  of  two  queues  in  tandem  where  the 
first  server  has  deterministic  service  times,  the  second  server 
has  general  service  distribution,  and  the  arrival  process  is  an 
arbitrary  renewal  process.  The  result  is  extended  to  the  case  of 
n queues  in  tandem  where  all  the  servers  except  the  last  one 
have  constant  service  times. 

^sing  a definition  of  variability  of  random  variables,  it  is 
proven  that  for  a tandem  queueing  system  with  n stations  in 
series,  where  each  station  can  have  either  one  server  with  an 
arbitrary  service  distribution  or  a number  of  constant  servers  in 
parallel,  the  expected  total  waiting  time  in  system  of  every  cus- 
tomer decreases  as  the  variability  of  interarrival  and  service 
distributions  decreases.  A new  sufficient  condition  for  customer 
average  delay  to  be  smaller  (larger)  than  time  average  delay  in 
single  server  queues  is  also  given 

A heterogeneous  arrival  single  server  queueing  loss  model  is  also 
analyzed,  where  the  arrival  process  is  a nonstationary  Poisson 
process  with  an  intensity  function  whose  evolution  is  governed  by 
a two-state  continuous  time  Markov  chain.  The  explicit  loss  formu- 
la for  this  model  is  obtained.  In  a special  case,  it  is  shown 
that  as  the  arrival  process  becomes  more  stationary  the  loss  de- 
creases. For  single  server  loss  systems  with  renewal  arrivals, 
counterexamples  are  given  to  show  that  regularity  of  arrival  and 
service  distributions  do  not  work  to  good  effect  in  general.  Two 
sufficient  conditions  for  it  to  be  true  are  given. 
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CHATTER  1 


INTRODUCTION 

This  thesis  is  concerned  with  two  types  of  problems:  bounds  for 
expected  delays  in  some  tandem  queueing  systems  which  are  very  difficult 
to  analyze  exactly  in  general  and  the  effect  of  regularity  of  arrival 
process  and  service  times  on  the  performance  of  some  queueing  systems. 

A tandem  queue  is  a number  of  service  facilities  in  series. 
Customers  arrive  according  to  a renewal  process.  Upon  arrival,  each 
customer  goes  to  the  first  station  and  requires  a random  amount  of 
service  time.  If  there  were  other  customers  present  at  the  time  of 
his  arrival,  he  joins  the  end  of  the  queue  and  waits  for  service. 

The  order  of  service  of  customers  is  first  come  first  served.  After 
being  served  at  the  first  station,  he  goes  to  the  second  station. 

Each  station  operates  in  a similar  fashion  but  may  provide  different 
type  of  service  in  general.  Every  customer  has  to  go  through  all 
stations  according  to  a prefixed  order  and  leaves  the  system  after 
finishing  service  at  the  last  station. 

Most  known  results  in  queueing  theory  are  for  queues  with  renewal 
arrivals.  The  fact  that,  in  general,  the  output  process  of  a queueing 
station  is  not  a renewal  process  makes  tandem  queues  very  difficult  to 
analyze. 

A number  of  authors  have  tried  to  analyze  the  output  process  of 


I 


a queueing  system  with  the  hope  that  it  might  be  helpful  in  analyzing 
the  behavior  of  a subsequent  station. 
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Burke  [5],  [6]  showed  that  Che  output  of  a stationary  M/M/ s 
queue  is  a Poisson  process  with  the  same  rate  as  the  arrival  process. 
Hence  In  a tandem  queue  with  Poisson  arrival  and  exponential  service 
times,  the  waiting  times  at  each  station  can  be  computed  by  solving 
a standard  M/M/1  queue.  Furthermore,  it  has  been  shown  by  Reich  [16], 
[17]  that  the  waiting  times  of  a customer  at  different  stations  are 
independent  of  each  other. 

For  a tandem  queueing  system  where  all  the  servers  have  constant 
service  times  except  maybe  one  whose  service  times  are  random  and 
larger  than  that  of  others  with  probability  one,  Friedman  [8]  showed 
that,  for  any  arbitrary  arrival  process,  the  epoch  at  which  each 
customer  departs  the  system  16  independent  of  the  order  of  servers. 

If  we  are  allowed  to  choose  the  order  of  servers  in  tandem  queues 
where  the  service  times  for  different  servers  are  nonoverlapping,  Tembe 
and  Wolff  [23],  [24]  showed  that  the  total  waiting  time  in  system  is 
stochastically  minimal  when  the  servers  are  in  decreasing  order  of 
their  service  times.  They  also  proved  that,  for  two  queues  In  tandem 
with  one  constant  server  and  one  arbitrary  server,  it  is  better  to  put 
the  constant  server  in  the  first  station  for  any  arbitrary  arrival 
process. 

Except  in  these  instances,  few  useful  results  are  known  concerning 
the  waiting  times  in  tandem  queueing  systems. 

In  view  of  the  difficulties  with  analyzing  tandem  queues,  it  is 


natural  to  look  for  some  approximations  and  bounds  for  various  quantities 
of  interest  in  these  models.  In  Chapter  2 and  3,  some  results  in  this 
direction  are  presented. 


Tembe  and  Wolff  [24]  obtained  an  upper  bound  for  the  expected 
waiting  time  for  a system  of  two  queues  In  tandem  where  the  first 
server  has  constant  service  times,  the  second  server  has  exponentlsl 
service  times,  and  the  arrival  process  to  the  system  is  Poisson. 

In  Chapter  2,  we  study  a more  general  system  with  a deterministic 
first  server  but  the  second  server  has  an  arbitrary  service  time  dis- 
tribution and  the  arrival  process  is  a general  renewal  process. 
Motivated  by  Tembe  and  Wolff '6  bound,  a conjecture  is  made  to  the 
effect  that,  for  a fixed  arrival  process,  the  expected  delay  In  front 
of  the  second  server  is  a decreasing  function  of  the  magnitude  of  the 
constant  service  time  at  the  first  station. 

If  the  conjecture  is  true,  then,  letting  the  service  time  at 
the  first  station  go  to  zero,  we  have  that  the  expected  delay  in  front 
of  the  second  server  is  bounded  above  by  the  corresponding  one  for  a 
system  without  the  first  server  or  just  a conventional  GI/G/1  queue. 

Well-known  upper  bounds  for  the  expected  delay  in  GI/G/1  queue 
are  readily  available.  Hence  the  conjecture  will  give  upper  bound  for 
the  expected  delay  in  front  of  the  second  server  in  the  above  system. 

We  have  been  unable  to  verify  the  conjecture  at  this  time.  How- 
ever the  upper  bound  which  can  be  obtained  in  this  fashion  is  shown 
to  be  valid . 

The  above  result  can  be  generalized  to  find  bounds  for  expected 
delays  in  tandem  queues  with  n stations  (n  >_  2)  in  series  where 
all  servers  except  the  last  server  whose  service  distribution  may  be 


arbitrary  have  constant  service  times  and  the  arrival  process  is 
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As  a by-product  of  the  analysis,  we  also  give  some  simple  bounds 
for  expected  Idle  periods,  variance  of  Interoutput  times,  and  some 
correlation  relations  among  various  quantities  In  some  tandem  queueing 
systems. 

In  Chapter  3,  ve  investigate  the  effect  of  regularity  of  the 
arrival  process  and  service  times  on  the  waiting  time  of  a customer 
in  some  tandem  queueing  systems. 

The  following  "conventional  wisdom"  16  generally  believed  to  be 
true  for  many  queueing  systems:  The  more  regular  the  arrival  process 
and  the  service  times  are,  the  better  the  system  performance  will  be. 

Wolff  [23]  gave  a summary  of  known  results  In  the  literature 
which  are  in  support  of  the  conventional  wisdom  and  also  presented  a 
number  of  systems  which  exhibit  contrary  behaviors. 

Stoyan  [18],  [21]  proved  a theorem  on  the  comparison  of  delays  in 
GI/G/1  queues.  He  gave  a definition  of  variability  which  is  very  useful 
in  this  context,  and  proved  that,  in  a GI/G/1  queue,  as  the  interarrival 
times  and  service  times  become  more  regular,  the  delay  decreases.  In 
this  thesis,  when  we  say  "more  regular,"  it  will  mean  more  regular  in 
Stoyan 's  sense  in  most  cases. 

His  result  can  readily  be  applied  to  give  a new  sufficient  condition 
for  customer  average  delay  to  be  smaller  (or  larger)  than  time  average 
delay  in  a GI/G/1  queue.  Improvements  of  some  known  bounds  by  Marshall 
[13]  for  the  expected  delay  in  certain  classes  of  GI/G/1  queues  also 
follow  easily. 

The  main  results  in  this  chapter  are  extensions  of  Stoyan 's  theorem 
to  some  tandem  queueing  systems.  It  is  shown  that  for  a tandem  system 
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with  n stations  in  series,  where  each  station  can  have  either  one 
server  with  an  arbitrary  service  distribution  or  a number  of  constant 
servers  in  parallel,  the  expected  total  waiting  time  of  every  customer 
decreases  as  the  interarrival  and  service  times  become  more  regular. 

The  results  are  quite  useful  for  bounding  purposes  since  we  can 
bound  the  expected  total  waiting  time  for  systems  which  are  difficult 
to  analyze  by  the  corresponding  quantity  for  easier  ones.  Two  examples 
of  this  sort  are  given. 

Under  stronger  assumptions  on  the  arrival  process  and  service 
times,  stronger  comparison  results  can  also  be  obtained. 

Specifically,  consider  a sequence  of  two  queues  in  tandem.  For 
two  such  systems  where  the  servers  at  the  first  station  have  the  same 
distribution  for  both  systems,  it  is  proven  that  the  system  with 
stochastically  larger  interarrival  times  and  stochastically  smaller 
(more  regular)  service  times  at  the  second  station  has  stochastically 
smaller  (more  regular)  delay  for  every  customer  in  front  of  the  second 
station. 

In  Chapter  4,  the  effects  of  regularity  on  some  single  server 
queueing  loss  models  are  examined. 

Fond  and  Ross  [22]  considered  a single  server  exponential  queueing 
loss  model  where  the  arrival  and  service  rates  alternate  between  two 
phases  (X^.u^)  and  (X2»p2^  • The  a®ounts  of  time  the  system  spends 
in  each  phase  have  exponential  distributions  with  rate  ca^  and  co^ 
respectively.  All  arrivals  who  find  the  system  busy  are  lost.  By 
solving  balance  equations,  they  showed  that  the  proportion  of  customers 
lost  is  a decreasing  function  of  c , i.e.,  the  more  stationary  the 
arrival  process  is,  the  smaller  the  proportion  of  customers  lost  will  be. 
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CHAPTER  2 

BOUNDS  FOR  SOME  TANDEM  QUEUEINC  SYSTEMS 


2.1  Definitions  and  Notations 

The  first  system  we  will  consider  is  a system  of  two  queues 
in  tandem.  We  will  denote  such  a system  by  GI/C^/1  ■*  Gj/l  where 
G^  , i « 1,2  is  the  distribution  of  the  service  times  at  station  i 
and  GI  means  that  the  arrival  process  is  a general  renewal  process. 
Similar  notations  will  be  used  throughout,  e.g.,  GI/D/1  ■*  G/l  , 

M/D/1  -*■  M/1  G/l  where  D , M represent  constant  and  exponential 
service  or  interarrival  times  respectively. 

By  "delay"  we  shall  always  mean  the  amount  of  time  a customer 
spends  waiting  in  queue  in  front  of  a service  station.  The  waiting 
time  of  a customer  is  equal  to  his  delay  plus  service  time. 

Unless  otherwise  specified,  it  is  assumed  that  the  queueing  process 
starts  at  time  zero  with  all  stati< ns  empty.  Hence  the  delay  of  the 
first  customer  will  be  zero. 

Let  the  subscript  n (e.g.,  Dfi)  refer  to  the  nCh  customer. 

The  following  notations  will  be  used: 


T = the  interarrival  time  between  the  nC^  and  (n  + l)t'1  arrival, 
n 


E(T  ) «=  ~ , Var  T 
n A n 


T 


Tn  • " 


1,2,  ...  , are  l.i.d. 


random  variables. 


S = service  time  of  the  nC^  customer  at  the  first  station, 
n 


E(S  ) « — , Var  S 
n n 


°S  < “ 


S , n * 1,2,  . . . , are  i.i.d. 
n 


random  variables. 
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R - service  time  of  the  nth  customer  at  the  second  station, 
n 

1 2 

E(R  ) “ ~ * Var  R ■ Or,  < “ • S , n - 1,2,  ...  , are  i.i.d. 
« 1*2  n k n 

random  variables. 

■ delay  in  queue  of  the  ntl1  customer  in  front  of  the  first 
station . 

* th 

D “ delay  in  queue  of  the  n customer  in  front  of  the  second 
n 

station. 

P1  “ *^wl  » x < • 

p2  “ */v2  » ^ < V*2  • 


W = D + S . 

n n n 

* * 

W - D + R . 

n n n 


2.2  A Conjecture 

Tembe  and  Wolff  [24]  proved,  among  other  things,  a useful  result 
concerning  the  optimal  order  of  servers  in  tandem  queues:  for  a tandem 
queue  with  two  servers  in  series,  if  one  of  the  server  has  constant 
service  time,  the  total  waiting  time  in  system  for  every  customer  is 
stochastically  smaller  when  the  constant  server  is  in  the  first  station. 
The  conclusion  is  true  for  any  arbitrary  arrival  process. 

By  reversing  the  order  of  servers,  the  above  was  used  to  find  an 
upper  bound  for  the  expected  stationary  waiting  time,  E(W*)  , in  front 
of  the  second  server  in  a M/D/1  -*■  M/1  system.  Specifically,  they 
obtained 


E(W*) 


(2.2.1) 
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Notice  that  the  right-hand  side  is  just  the  expected  waiting  time  of 
a M/M/1  queue.  The  intuitive  explanation  of  this  is  that  the  constant 
server  at  the  first  station  tends  to  make  the  arrival  process  to  the 
second  station  more  regular  than  it  would  be  without  the  first  server. 

It  is  felt  that  the  validity  of  the  above  explanation  should  not 
depend  on  the  distributions  of  the  arrival  process  and  service  time  at 
the  second  station.  Hence,  we  make  the  following 


Conjecture: 

In  a GI/D/1  -*■  G/l  system,  the  expected  delay  in  front  of  the  sec- 
ond server  is  smaller  than  it  would  be  if  there  were  no  first  server  at 
all. 

In  fact  we  believe  a stronger  conjecture  should  also  be  true, 
namely,  the  larger  the  constant  service  time  at  the  first  station,  the 
smaller  the  expected  delay  in  front  of  the  second  server  will  be. 

Thus,  the  first  conjecture  follows  by  letting  the  service  time  at 
the  first  station  go  to  zero. 

If  the  conjecture  were  true,  then  It  follows  from  well-known 
upper  bound  for  the  expected  delay  in  GI/G/1  queues  by  Kingman  [9], 

[10]  and  Marshall  [13]  that 


(2.2.2) 


where  denotes  the  stationary  delay  of  a customer  in  front  of  the 

second  station. 


We  have  been  unable  to  verify  the  conjecture  at  this  time.  How- 
ever, it  is  proven  in  the  next  section  that  (2.2.2)  is  indeed  valid. 
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2.3  Bounds  for  the  Expected  Delay  in  Front  of  the  Second  Server 


It  is  well-known  that 


(2.3.1) 


Vl  ' ■“  |0>D„  + Sn  - V 


or  equivalently, 


(2.3.2) 


D . - X - D + S -T 
n+1  n n n n 


where  X ■ -min  [0,D  + S - T ] . 

n n n n 

We  will  first  give  relations  similar  to  (2.3.1)  and  (2.3.2) 
and  then  use  them  to  prove  (2.2.2)  by  an  approach  parallel  to  Kingman 
and  Marshall's. 

To  start  with,  we  do  not  assume  that  the  service  times  of  the 

first  server  are  deterministic.  Let  e^  be  the  epoch  at  which  the  nth 

* 

customer  arrives.  Observe  that  e+D+S+D+R  is  the  epoch  he 

n n n n n 

leaves  the  system  and  e + T + D + S , , is  the  epoch  the  (n  + l)th 

n n n+i  n+i 

customer  departs  the  first  station.  Thus 


(d  +S  +D*  + R -T  - D . - S , 
I n n n n n n+1  n+1 

|o  otherwise. 


, if  D + S + D + R > T + D . , + S , , 
n n n n n n+1  n+1 


Equivalently, 


(2.3.3)  D , « max 
n+l 


|o,D  + S + D*  + R -T  - 
L n n n n n 


Dn+1  " Sn+l]  * 


Let  Y “ -min  (0,D  + S +D  + R -T  - D..  - S.,"l,we  have 

n L n n n n n n+1  n+lj  ’ 
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(2.3.4)  D . - Y - D + S + D +R  - T - D . - S . . 

n+1  n n n n n n n+1  n+1 


In  the  derivations  below,  we  will  assume  stationarity  of  the 
processes  {D^.n  ±_  1)  and  {Dn,n  — * 

Taking  expectations  in  (2.3.4),  we  have 


(2.3.5) 


E(Y  ) ■ E(T  ) - E(R  ) 
n n n 


f (1  - P2>  * 


Note  that  (2.3.5)  can  be  used  to  find  a lower  bound  for  the 
expected  idle  period,  E^)  , at  the  second  station  by  using 
E(Y^)  * SqE^)  <_  E^)  where  a^  = the  probability  that  a stationary 
customer  departing  from  the  first  station  finds  the  second  station 
empty. 

Let  B *=R  -T  - S , (2.3.4)  can  be  written  as 
n n n n+1 

(d  ..  + D*  \ - Y - (D  +D*)+S  + B . 

\ n+1  n+1 ) n\n  n J n n 

★ 

Squaring  both  sides  and  noting  that  D^+^Y^  = 0 gives 

(d  + D*  )2  - 2Y  D + Y2  - (d  + D*\2  + S2  + B2 
\ n+1  n+1  j n n+1  n \ n nf  n n 

+ 2|D  + D*)S  + 2|D  + D*  )B  + 2S  B . 

\n  n/n  \n  n/n  nn 


Note  that  the  pairs  D and  S , D and  B , D and  B , 

n n n n n n 

S and  B are  independent.  After  taking  expectations  and  can- 
n n 


celling,  we  get 
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(2.3.6) 


e(y2^  - 2E(Y  D ..)  - E (s2  ) + E(B2)  + 2E(D  )E(S  ) 
\ n I n n+i.  \n/  \ n / n n 

+ 2E|D*S  ) + 2E(D  )E(B  ) + 2E(S  )E(B  ) . 

\ n n/  n n n n 


Substituting  E(YnDn+1)  - E(Yn)E(Dn+1)  + Cov  [Yn»Dn+l]  , 


E^SnDnj  » E (SR)E^Dn  j + Cov  £sn*DnJ  an<i  n°ting  that 

E(S  ) + E(B  ) - -E(Y  ) , (2.3.6)  becomes 
n n n 


n n+1 


'{‘I)  - 2 C-  [VVl]  - E(Sn)  + E(BJ)  - 2E(D»)E<V 

+ 2 Cov  [s  ,D* 1 + 2E(S  )E(B  ) . 

L n n J n n 


After  some  rearrangements,  we  have  proven  the  following 


Theorem  2.3.1: 


Under  stationary  conditions, 


(2.3.7) 


-5-7^ ^ < 2o2  + o2  + o2  - Var  Y 

2(1  - p2)  | S R T e 


+ 2 Cov  [Y  ,D  . , ] + 2 Cov 
n n+1 


M| 


Remark: 


In  the  derivations  of  (2.3.5)  and  (2.3.6),  we  needed  the  condi- 


K) 


tions  that  EID  1 < « and  E|D  1 < «■  in  order  to  make  cancellations. 


It  is  assumed  that  these  cancellations  are  valid.  In  practice,  this 

assumption  is  usually  justified. 

The  unknowns  in  (2.3.7)  are  Var  Y , Cov  [Y  ,D  , and 

n n n-rj. 

Cov  ^n,DnJ  . However,  we  always  have  Var  Y^  >_  0 and,  in  the  case 


t] 
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of  constant  first  server,  Cov  = 0 . Hence,  we  only  have 

to  find  an  upper  bound  for  Cov  tYn»Dn+^J  • In  order  to  do  that, 
the  concept  of  association  of  random  variables  is  used.  It's 
definition  and  basic  properties  are  given  in  Esary,  Proschan  and 
Walkup  [7].  For  convenience,  a summary  is  given  here: 


Definition: 

Random  variables  V V are  "associated"  if 

1 n 

Cov  [T(V)  ,MV)]  >_  0 for  all  pairs  of  increasing  binary  functions 
r , A . We  shall  also  say  the  vector  V ■ (V, , . . . , V ) is 
associated. 

Associated  random  variables  have  the  following  properties: 


(PI):  Any  subset  of  associated  random  variables  are  associated. 
(P2):  The  set  consisting  of  a single  random  variable  i6 
associated. 

(P3):  Increasing  functions  of  associated  random  variables  are 
associated. 

(P4):  If  two  sets  of  associated  random  variables  are  independent 
of  each  other,  then  their  union  is  a set  of  associated 
random  variables. 

(P5):  Independent  random  variables  are  associated. 

(Note:  this  follows  immediately  from  P2  and  P4). 

(P6):  If  random  variables  V.,  ...,  V are  associated,  then 

J.  n 

Cov  [f(V),g(V)l  0 for  all  increasing  functions  f 

and  g . 

Before  proving  the  next  theorem,  we  need  the  following 
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D , , - max 
n+1 


[°.  „i„  [d„  + T.  ' So+1,D*  + *„  - Sn+1]] 


everywhere . 

The  proofs  of  (2.3.9)  and  (2.3.10)  are  similar.! 

Unless  otherwise  stated,  we  will  assume  from  here  on  in  this 

section  that  S = s , Vn  > 1 , where  s is  a constant,  i.e.,  the 
n 

first  server  has  constant  service  time  s . In  this  case,  (2.3.8) 
and  (2.3.9)  reduce  to 


(2.3.11)  Dn+1  - max 


(2 


|0,  min  [d  + D*  + R - T ,D*  + R - si  1 
L’  [ n n n nn  n JJ 

.3.12)  Y = -min  fo.D^  + D*  + R - T D*  + R - s]  . 

n [_nnnnnnj 


We  proceed  to  the  following  important 


Theorem  2.3.3: 


For  all  n > 1 , D and  D are  associated  random  variables. 
— n n 


Proof : 


The  proof  is  by  induction. 

* * 

Obviously,  and  are  associated  since  ■ 0 . 

Assume  D and  D are  associated.  It  is  clear  from  (2.3.1) 
n n 

* 

and  (2.3.11)  that  both  Dn+1  and  Dn+1  are  increasing  functions 

of  the  vector  |D  ,D*,-T  ,R  ) . Now,  -T  and  R are  independent 
\ n n n n/  n n 

of  each  other  and  ID  ,D  | . Therefore,  it  follows  by  induction 

' 11  / ' * 

hypothesis  and  P4  that  |Dn,Dn,-Tn,R  \ is  associated.  Hence,  by  P3, 

* 

D and  D are  associated.! 

n+1  n+1 
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Lemma  2.3.4: 


For  all  n > 1 , -Y  and  D are  associated  random  variables. 
— n n+1 


Proof : 


It  follows  from  (2.3.1)  and  (2.3.12)  that  both  D , and  -Y 

n+1  n 

are  increasing  functions  of  the  vector  |Dn,Dn*~^n,^n j • But,  by 

Theorem  2.3.3,  P4,  and  P5,  [D  ,D  ,-T  ,R  ) is  associated.  Hence, 

’ \ n’  n’  n’  n/ 


the  conclusion  follows  from  P3.I 


Remark : 


As  n goes  to  infinity,  -Y  and  D ,,  converge  in  distribution 

n n+l 

to  their  stationary  distributions  respectively.  Since  (-Y^D^)  *8 

associated  for  all  n 1 , it  follows  that  the  stationary  distribu- 
tions of  -Y  and  D are  also  associated, 
n n 


We  are  now  ready  for  the  following  main 


Theorem  2.3.3: 


For  a stationary  GI/D/1  •+•  G/l  system, 


/ *\ 

: D < 
\ n/  — 


\°T  + °r) 
2(1  - P2) 


where  equality  holds  for  D/D/1  D/1  systems. 


Proof : 


If  S = s , Vn  > 1 , then  o_  * 0 and  Cov 
n — s 

Hence,  (2.3.7)  specializes  to 


[V>n] 
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\°l  + °R  - W"  \ * 2 C°*  IV ’Vl1} 

|o^  +o^+2  Cov  IYn»Dn+^]|  since  Var  Yr  0 


Cov  [Y  ,D  , ] <0  by  Lemma  2.3.4.B 
n n+1  — 


Remark: 

Bounds  for  the  variance  of  the  interdeparture  time  from  the 

second  station  can  also  be  obtained.  Let  t * time  between  the  nth 

n 

and  (n  + 1)  departure  from  the  second  station.  It  follows  from  the 

fact  that  t = Y + R that  E(t  ) - 1/X  . Since  Y and  R 

n n n+1  n n n+1 

are  independent,  we  have 


(2.3.13) 


Var  t = Var  Y + Var  R 

n n n+1 


— ^R 


In  the  case  = s , Vn  >_  1 , (2.3.7)  becomes: 


Var  Y = a2  + ol  + 2 Cov  [Y  ,D  . . ] - 2e(d*W(Y 
n R T n n+1  \ n/  n 


) 


(2.3.14) 


< al  + a2  + 2 Cov  [Y  ,D  since  E(D*\e(Y  ) > 

— R T n n+1  \ n)  n — 

2 2 

< o_  + o since  Cov  [Y  ,D  ..]  < 0 . 

— K i n n+1  — 


Combining  (2.3.13)  and  (2.3.14)  gives 


°R  £ Var  Tn  < 2c2r  + c2  , 


where  equality  holds  both  sides  for  D/D/1  •+  D/1  systems. 


j 
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2.4  A Generalization 


Consider  the  following  system  of  n(n  >_  2)  queues  in  tandem. 
Customers  arrive  according  to  a general  renewal  process  where  the  inter- 
arrival times  T , n * 1,2,  ...  are  i.i.d.  random  variables  and  have 
n 

1 2 

an  arbitrary  distribution  with  E(T  ) *=  — , Var  T « o_  . The  first 

n a n T 

through  (n  - l)1"*1  stations  have  constant  servers  with  deterministic 
service  times  , j = 1,2,  ...,  n-1  , and  the  nth  server  has  an 

2 

arbitrary  service  distribution  H with  E(H  ) « 1/y  , Var  H « o„  . 

n n n n h 

Customers  are  served  according  to  first  in  first  out  rule  at  each 
station.  Under  stationary  conditions,  an  upper  bound  for  the  expected 
total  waiting  time  in  system  of  a customer  is  found  below. 

Friedman  [8]  has  shown  that  if  all  the  servers  in  a tandem  queueing 
system  have  constant  service  times,  then,  for  any  arbitrary  arrival 
process,  the  epoch  at  which  every  customer  leaves  the  entire  system 
does  not  change  with  changes  in  order  of  servers.  Now,  if  we  consider 
the  first  through  (n  - 1)^  server  as  the  first  subsystem  and  the  nC^ 
server  as  the  second  subsystem,  then  our  system  becomes  two  subsystems 
in  tandem. 

In  the  first  subsystem,  since  the  total  waiting  time  is  independent 
of  the  order  of  servers,  we  can  rearrange  them  such  that  their  service 
times  are  in  decreasing  order.  Let  d^j  , j = 1,  . . . , n - 1 , be  the 
delay  in  front  of  the  station  after  the  rearrangement,  then 
d.  i = 0 for  all  j = 2,  ...,  n-1  , and,  by  Kingman  and  Marshall's 


bound , 


HiiHj 


E(d[i])  -Tonh y where 

Vj  * 1,  ....  n-1}  . 


k c (i  | i e {1, 


n-1}  , 
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Let  (t^  , 1 “ 1,2,  ...}  be  any  realization  of  the  departure 
epochs  from  the  first  station  after  rearrangement,  then  it  16  easy 
to  see  that  the  departure  epochs  from  the  first  subsystem  are 


| "j1  | 

'i+  I “j  • 


r a-1 


Now,  if  we  compare  our  system  with  a GI/H^/1  -*■  G/l  system  (i.e. 


>th 


delete  the  second  through  (n  - 1)  stations),  the  above  observation  says 


that  the  arrival  processes  to  the  last  station  in  the  two  systems  differ 
n-1 


by  a constant  £ H - H,  . Hence  every  customer  will  experience  the 

j-1  J 

same  amount  of  delay  in  front  of  the  last  station  for  any  realization  of 
the  processes  in  both  systems.  Therefore,  by  Theorem  2.3.5,  the  expected 


2 ^ 2 ' 
M°t  °u 

delay  in  front  of  the  last  station  is  bounded  above  by  n — s 

2(1  - X/u  ) 


We  have  proven  the  following 


Theorem  2.4.1: 


*°T  n_1  »(°T  + °h)  1 

E(w)  1 2U  - + lml  -j  * * rn 


where  W = total  waiting  time  in  system  of  a stationary  customer. 


Again,  equality  holds  if  both  T and  H are  deterministic. 

n n 


Remark: 


(i)  Friedman's  result  holds  even  if  the  deterministic  servers 
have  multiple  parallel  channels.  Theorem  2.4.1  can  be 
easily  modified  to  this  case. 

(ii)  Tembe  and  Wolff's  bound,  (2.2.1),  can  also  be  generalized 
in  exactly  the  same  way. 
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2.5  Some  Correlation  Relations 


In  this  section,  we  Investigate  some  properties  of  stationary 
covariances  and  their  relationships.  We  begin  with  the  following 
Intuitively  obvious  theorem: 

Theorem  2,5.1: 

In  GI/G/1  G/l  systems,  Cov  Is  ,D  1 < 0 and  Cov  [S  , ,Y  ] > 0 

L n n J — n+l  n — 

for  all  n ^ 1 . 

Proof : 

* 

Conditioning  on  D , D , , R , , S , , T , , we  have 
n-i  n-i  n— i n— l n— ± 


+ Cov 


Cov  [s  ,D*1  - E < Cov  [s  ,D*  | D . ,D*  ,R  . ,S  , ,T  .1/ 
u nj  | L n n ' n-1  n-1  n-1  n-1  n-1  Jl 

HS  | D . ,D*  . ,R  ,,S  . ,T  ,\,e(d*  I D ,.R  ,.S  . ,T  .\\ 

n 1 n-1  n-1  n-1  n-1  n-1 )’  \ n 1 n-1  n-1’  n-1  n-1  n-l/i 


Since  S is  independent  of  D , D , , R , S , T , 
n n-i  n-i  n— l n-i  n-i 

It  follows  e/s  Id  , ,D  , ,R  , ,S  , ,T  - E (S  ) , a constant. 

^ n 1 n-1  n-1  n-1  n-1  n-lJ  n 

Thus  the  second  term  above  is  zero. 

Cov  fs  ,D  | D , ,D  . ,R  , ,S  , ,T  ,1  < 0 follows  from  the 
L n n 1 n-1  n-1’  n-1  n-1’  n-lj  — 

* * 
fact  that  D is  a decreasing  function  of  S given  D , , D , , 
n n n-1  n-1 

R , , S , T (see  (2.3.8)).  Therefore,  Cov  |S  ,D  1 < 0 . 
n-l  n-l  n-i  L n nJ  — 

Using  a similar  argument,  Cov  [S  . ,Y  ] > 0 follows  by 

n+l  n 

* 

conditioning  onD  , D , S ,R  , T . ■ 
n n n n n 

* 

In  the  case  of  deterministic  first  server,  D and  D are 

n n 


associated  for  all  n > 1 (Theorem  2.3.3).  We  have 


Theorem  2.5.2: 


Each  of  the  following  vector  is  associated  for  all  n ^ 1 : 

<-x„-D„>  * (-VV)  • <D»-Vi>  • (-VDn)  • 'VV  • (-VVi)  • 

K-Cl)  • (Dn,Dn+l)  ' (\'Vl)  • a"d 


Proof : 


Similar  to  the  proof  of  Lemma  2. 3. 4. 1 


Next,  we  will  derive  some  relationships  and  bounds  for  various 
stationary  covariances. 

Observe  that  (2.3.4)  can  also  be  written  as: 


(2 


.5.1)  (d  ..  + S ..  + D*  ) - Y * (d  + S + D*)+R  -T  . 
\ n+1  n+1  n+1  f n\n  n n J n n 


Now,  if  we  apply  the  same  procedures  to  (2.3.10)  and  (2.5.1)  as 
in  the  derivation  of  (2.3.7),  we  can  get  (omitting  the  details): 


(2.5.2)  E 


I *\  1 

\ n ) 2E(Y  ) 


Var  X +o^  + Op  - Var  Y -2  Cov 
I n R S n 


[Xf<]J 


and 


(2-5-3)  E(D»)-s^y  K+0?-v“  Y„+2  Cov  <WiI  + 2 Cov  IY„-Vi 

Therefore,  by  comparing  (2.3.7),  (2.5.2),  and  (2.5.3),  we  obtain 

Theorem  2.5.3: 

For  a stationary  GI/G/1  -*•  G/l  system. 
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(2. 5. A)  Var  X -2  Cov  [x  ,D*|-o^  + o^+2  Cov  [Y  ,D  + 2 Cov  fs  ,D*1  , 
n L n nj  a i n n+i  [_  n nj 

(2.5.5)  „l  ♦ Cov  [sn.D*J  - Cov  IT^I  . 

(2.5.6)  Var  °S  ' 2 Cov  [x„,I>*]  - 0*  + 2 Cov  [Y^D^l  + 2 Cov  [V^S^l  . 

Corollary  2.5.4; 

Under  stationary  conditions, 

-o2  <Cov  [s„.»;]iO 

0 i C0V  1 WlJ  i »S 

for  GI/G/1  -*•  G/l  systems.  Equality  holds  both  sides  if  S 's  are 

n 

deterministic . 

Proof : 

By  Theorem  2.5.1  and  (2.5.5).B 

Corollary  2.5.5: 

Under  stationary  conditions, 

Var  Xn  - 2 Cov  [x„.D*]  -0^+2  Cov 

if  S 's  are  deterministic, 
n 

Proof : 

By  (2.5.4)  and  (2.5.6).B 


J 
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Corollary  2.5.6: 

If  the  arrival  process  is  Poisson  and  S » s Vn  > 1 , 

n — 

then,  under  stationary  conditions 

- T i c”  [vC]  i 0 
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Proof : 


It  is  well-known  that  W and  W are  independent  (see 

n n 

Reich  [17]).  Since  X - -min  [0,W  - T ] and  T is  independent 

n n n n 

* * 
of  W , it  follows  from  the  above  fact  that  X and  W are 
n n n 

* 

independent.  This,  in  turn,  implies  X and  D are  independent. 


Therefore,  Cov 


[vl] 


n 


n 


0 . 


Noting  that  Var  X ■ , (2.5.4),  (2.5.6)  simplify  to 

X Wl 

(2.5.7)  and  (2.5.8)  respectively.® 


The  exact  value  of  Cov  [Y  ,D  ,,]  can  be  found  by  using  the 

n n+l  ° 

* 

fact  that  and  are  independent  for  the  M/M/1  ■+■  M/1  system. 

Explicitly,  we  have 


E<Y„Vl>  * E(Y„Vl  I Vl  - °>P<Vl  * °> 

■ E(Y„  I Vl  * °>E‘Vi  I Vl > 0)p{Vi > 01  • 


where  we  have  used  the  assertion  that  (Y 


D >0)  and 
n+l 


(Dn+^  | > 0)  are  independent.  The  assertion  follows  since 


<Y„  I Vi » °> 


max 


[°'Vi  - <]  • <Vi  I Vl - 0> 


(W  - T 
n n 


W > T ) , and  the  vectors  |S  ,W  \ and  (W  ,T  ) are 
n n \ n+l  n / n n 


independent  of  each  other.  Now, 


E(Vi  I Vl > 0,P1Vi > 01  ■ E(Vi> 


Pl^l  - x)  ’ 


I 


E(Vn  | Dn+1  > 0) 


( r 

*1 1 

(max 

s . . - w s > 

n+1  n n+1 


-E( 

1 ( -2  - > \ 
"l  W + "2  - V 


<Hs»+i  > <> 


Hence,  E(Y  D 

n n+3  y 


*1  (“l  * “2  ' * X>)  ’ 


and 


C0V  lYn*Dn+l] 


E(YnVl)  - E(Yn)E(Dn+l) 


(2.5.9) 


x(u2  - X) 


(li2  ~ 
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y^^  - xHw^  + y2  - X)  Xy2 


y1(y1  - x) 


(y2  - x)(w1  + p2) 


wly2^i  + w2  " 


Therefore,  Cov  fs  ,D  1 and  Cov  [Y  ,S  can  be  obtained  via 

L n nj  n n+1 

(2.5.9)  and  Corollary  2.5.7. 
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in  parallel  h,  or  a number  of  rri 
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We  wil 1 start  with  c 

some  preliialnar^es^ 
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aU  distribution 


if 


FOr  tV°  dlStr—  P M c 

G * we  say  that  p ? G 
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/f(x)dF(x)  < J f (x)dG(x) 
for  all  Increasing  convex  function  f . 

Definition: 


For  two  arbitrary  random  variables  X and  Y , we  say  that 
v 

X <_  Y if  their  corresponding  distribution  functions  satisfy  the 
above  ordering. 

A list  of  some  useful  properties  of  this  ordering  follows: 

(a)  If  X < Y , then  E(x)  < E(Y)  . 

(b)  F G if  and  only  if 


00  00 

J [1  - F(x)]dx  <_j  [1  - G(x)]dx  Vt  . 

t t 

(See  Theorem  14  of  Bessler  and  Veinott  [2].) 

v 

(c)  If  F and  G have  the  same  mean,  then  F G if  and  only 


/f(x)dF(x)  < /f(x)dG(x) 

for  all  convex  function  f (see  Ross  [19]). 

(d)  Let  F and  G be  two  distribution  functions  with  the  same 
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Before  giving  the  next  property,  we  need  the  following  definition: 
Definition: 

A nonnegative  random  variable  X with  finite  mean  is  said  to 
have  NBUE  (NWUE)  distribution  if 

E(X  - t | X > t)  E (X)  Vt  > 0 . 

(e)  Let  X be  a nonnegative  random  variable,  and  Y be  an 

exponential  random  variable  with  the  same  mean  as  X . 

Then  X < Y if  and  only  if  X has  NBUE  (NWUE)  dis- 

(1) 

tribution.  (See  Theorem  4.6  of  Marshall  and  Proschan  [15]). 

It  should  be  noted  that,  for  two  random  variables  X and  Y with 

v 

the  same  mean,  property  (c)  implies  that  Var  X <_Var  Y if  X ^ Y . 
Hence,  this  ordering  indeed  orders  the  relative  variabilities  of 
random  variables  with  the  same  mean. 

We  now  state  the  following  theorem  of  Stoyan: 

Theorem  3.2.1: 

Let  F^/G^/l  and  ¥^^2^  ^e  two  s^n8^e  server  queues  with  inter- 
arrival time  distributions  , i ■ 1,2  , and  service  time  distribu- 

00  00  ^ 

tions  G.  , i ■ 1,2  . If  / xdF. (x)  * / xdF, (x)  , F.  <_  F,  , and 

1 0 0 1 1 
v 1 v 2 i 

G^  G^  , then  <_  , where  D^  , i ■ 1,2  , are  the  stationary 

delays  of  systems  F^/G^/l  , i * 1,2  , respectively. 
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3.3  Customer  Average  and  Time  Average  Delays 

In  this  section,  some  relations  between  customer  average  and 

time  average  delays  for  some  classes  of  single  server  queues  are 

investigated.  By  "customer  average  delay,"  we  mean  the  average  amount 

of  work  in  system  as  seen  by  an  arriving  customer,  and  time  average 

delay  is  defined  as  the  average  over  time  of  the  amount  of  work  in 

system  at  an  arbitrary  time  point  t . It  should  be  noted  that  the 

amount  of  work  in  system  at  time  t is  equal  to  the  amount  of  time  a 

customer  has  to  wait  if  he  arrives  at  time  t . 

Marshall  and  Wolff  [14]  showed  that  the  time  average  delay  is 

greater  than  the  customer  average  delay  in  a single  server  queue  if 
2 

po  <1  where  p is  the  service  rate  and  o is  the  coefficient  of 
a — a 

variation  of  the  interarrival  time. 

In  the  next  theorem,  we  use  Stoyan's  theorem  to  obtain  a new 
sufficient  condition  for  time  average  delay  to  be  greater  (smaller) 
than  customer  average  delay  in  GI/G/1  queues. 

Let  E(V)  and  E(D)  be  the  time  average  delay  and  customer 
average  delays  in  a GI/G/1  queue  respectively.  Brumelle  [4]  showed  that 
the  following  relation  holds: 


E (V)  * pE  (D)  + 


AE(S2) 


where  X is  the  arrival  rate,  E(S)  and  E(S  ) are  the  first  and 
second  moments  of  the  service  time  distribution,  and  p «=  AE(S)  . It 
follows  that  E(V)  ^ 


E(D)  if  and  only  if 
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(3.3.1) 


E(D) 


<_  XE(S2) 

(>)  2(1  - p)  ' 


Now,  compare  the  original  GI/G/1  system  with  a M/G/l  system 
where  the  service  distribution  stays  the  same  but  the  arrival  process 
is  Poisson  with  rate  X . Notice  that  the  right-hand  side  of  (3.3.1) 
is  just  the  expected  delay  for  the  M/G/l  system.  Therefore,  the  next 
theorem  follows  from  property  (e)  of  variability  ordering  and  Theorem 
3.2.1. 

Theorem  3.3.1: 

In  GI/G/1  queues,  customer  average  delay  is  smaller  (greater)  than 
time  average  delay  if  the  interarrival  times  have  i.i.d.  NBUE  (NWUE) 
distributions . 

Denote  a single  server  queueing  system  with  i.i.d.  NBUE  (NWUE) 

interarrival  time  distributions  by  a NBUE/G/1  (NWUE/G/1)  system. 

As  a consequence  of  above,  Marshall's  bounds  [13]  for  the  expected 

2 

delays  in  NBUE/G/1  and  NWUE/G/1  queues  can  be  improved.  Let  ofl 
2 

and  a be  the  variances  of  the  interarrival  time  and  service  dis- 
g 

tributions  respectively.  We  have 


(i)  For  NBUE/G/1  queues  with  p < 1 , 


(3.3.2) 

(ii) 


■“  [t-J  - -Sr21]  i i ■*”  [J>  irf^y]  • 

For  NWUE/G/1  queues  with  p < 1 , 


max 


XE(S2)  I 

2(1  - p)J 


< E(D) 


< J 


(1  + P) 
2X 


(3.3.3) 


./  2 , 2 

Al  O + O ® 

where  J = 2(1 ) an(*  t'ie  solution  of  f [1  - K(u)]du  , 

-x 

x >_  0 , K is  the  distribution  of  the  difference  of  a service  time 
and  an  interarrival  time. 

Letting  C = Xa  , we  have 

3 3 


9 r2  n 2 

1E(S2)  ca  ~ P 

2(1  - p)  = 2A (1  - p) 


Thus , 


. j . Lii  „ c .1 

2(1  - p)  2\  a 


Therefore,  the  upper  and  lower  bounds  converge  to  each  other  as 
Ca  - 1 in  both  (3.3.2)  and  (3.3.3).  This  implies 


E(D)  - as  ca  - 1 • 


Hence  | E (D)  - E(V)  | -*•  0 as  1 for  both  NBUE/G/1  and 


NWUE/G/1  queus. 


Remark : 


Marshall's  bound  for  IFR/G/1  and  DFR/G/1  queues  can  also  be 
improved  in  the  same  way. 
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3.4  Comparison  of  Waiting  Times 

We  now  proceed  to  generalize  Theorem  3.2.1  to  tandem  queueing 
systems . 

Consider  a queueing  system,  , with  n stations  in  series. 

There  is  only  one  server  at  each  station.  Customers  arrive  according 
to  a renewal  process  with  interarrival  time  distribution  F . The 
service  times  at  station  j have  distribution  , j = 1,2,  ...,  n . 

Let 

T^  = interarrival  time  between  the  (k  - l)th  and  customer, 
k *=  1,2 

= service  time  of  the  k^  customer  at  station  j , 

j * 1,2,  . . . , n , k = 1,2 

Zj|  = epoch  at  which  the  k*"*1  customer  leaves  station  j , 

k 

j = 1,2,  ....  n , k = 1,2 Z°  = l T . 

i=l  1 

We  make  the  assumption  that  the  system  starts  operation  at  time 
zero  with  no  customer  in  system,  and  customers  are  served  in  the 
order  of  their  arrival  at  each  station. 

The  following  lemmas  are  essential  in  the  proofs  below: 

Lemma  3.4.1; 

Let  . Y , i ■ 1,2,  ....  n be  independent  random  variables, 
v 

Then,  for  all  i - 1,2,  ...,  n if  and  only  if 

f (X1,X2 Xn)  - f (Y1’Y2 V 


for  all  increasing  convex  function  f . 
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Proof : 

The  if  part  is  trivial,  so  we  only  have  to  prove  the  other 
direction.  The  proof  is  by  induction  on  n . 

The  lemma  obviously  holds  for  n « 1 by  definition  of  vari- 
ability. Suppose  it  is  true  for  k . Denote  F (F  | the 

k+1 \ Yk+l/ 

distribution  of  X^+1(Yk+^)  * Let  f and  8 be  arbitrary  increasing 
convex  functions.  Noting  that  increasing  convex  functions  of  increasing 
convex  functions  are  still  increasing  convex,  we  have 


E{g[f(X1,X2,  ...,  X^X^)]} 

E{g[f(X1,X2,  ...,  I “ Xk+l}dFXk+1(xk+l) 

- E{g[f(Y1,Y2,  ....  Yk»xk+1)i  I xk+1  “ Xk+l}dFXk+1*Xk+l) 
by  induction  hypothesis 

£ E{g[f(Y1,Y2,  ....  Yk>xk+1>]  I Yk+1  = xk+i}dFY  (xk+lJ 

k+1 

since  E{g[f(Y1,Y2,  ....  Yk>xk+1)]  I XR+1  = XR+1>  is 

V 

an  increasing  convex  function  of  and  — ^+1 

= E{g[f(Y1,Y2,  ....  Yk.Yk+1)]>  • ■ 


Lemma  3.4.2: 

Z^  is  an  increasing  convex  function  of  T^,  ....  T^,S^,  ....  S*. 
...,  S^,  ...,  for  all  j ■ 1,2,  ...,  n , k = 1,2,  ...  . 

Proof: 

Observe  that  max  ^Zj|  ^»zk  is  the  epoch  at  which  the  kth 
customer  enters  service  at  station  j . Hence 
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- max  [zj  1 ,Z^  J + for  all  J - 1,2,  . . . , n , k ■ 1,2,  ... 


O p 

Since  Z - 2 T , k « 1,2,  ...  , the  conclusion  follows  by 

k 1-1  1 

Induction  on  j and  k . ■ 


Now,  consider  a second  such  system,  , with  Interarrival  dis- 
tribution F and  service  distribution  at  station  J . Let 

all  wiggled  notations  (e.g.  Z^)  denote  the  corresponding  entitles 

for  system  . It  follows  from  Lemma  3.4.1  and  3.4.2  that 


(3.4.1) 


Z^  for  all  j - 1,2,  ....  n , k - 1,2,  ...  , 


v - 


if  F F and  for  all  J ■ 1 , 2 , ...,n. 

Let  Wk (Wk ) be  the  total  waiting  time  in  system  , we 

have 


Theorem  3.4.3: 


If  JxdF(x)  =y*xdF(x)  < °°  , F <_  F , and  G^  <_  G^  for  all 
j = 1,2,  ....  n , then  E(Wk)  1.  E(W^)  for  all  k - 1,2,  ...  . 


Proof : 


Wk  ' Zk  ' 2 Ti 

K K i-1  1 


w = z"  - y t. 

k k i^;1  i 


Taking  expectations  and  substracting , we  have: 
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E(V  - E<"k>  ■ e(z£)  - E(zl!)  i 0 • 

where  the  last  inequality  follows  from  (3.4.1).B 

If  the  arrival  processes  of  and  ^ are  the  same,  a 

stronger  relation  between  W and  W holds. 

Theorem  3.4.4: 

If  F(x)  = F(x)  for  all  x >_  0 , then 


Wk  <_  Wk  for  all  k - 1,2,  ...  . 


Proof : 

Let  f be  an  arbitrary  increasing  convex  function.  Condition- 

k 

ing  on  the  sequence  of  arrival  epochs  lt^}^=^  , we  have 

E)f<v  i i 

since 

Si  (Si)  ’ j = 1,2 n * 1 = 1,2 k * 

Unconditioning,  we  have  E[f(k'k)]  <_E[f(Wk)]  . ■ 

Remark: 

In  the  proof  of  Theorem  3.4.4,  we  only  needed  the  assumption  that 
k 

the  sequence  has  the  same  joint  distribution  for  both 


{t 


is  an  increasing  function  of 


systems.  Therefore,  the  interarrival  times  may  be  "dependent." 
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Another  interesting  question  one  might  ask  is:  Can  Theorem  3.2.1 
be  generalized  to  parallel  server  queues?  Ross  [19]  gave  a counter- 
example showing  that  similar  comparisons  do  not  hold  in  general. 
However,  we  have  the  following  theorem: 


Theorem  3.4.5: 


Let  A/S^/m  and  B/S2/m  be  two  parallel  server  queues  with  inter- 
arrival time  distribution  A and  B , deterministic  service  dis- 
tributions and  S£  , and  m servers  in  parallel  respectively. 

If  JxdA(x)  = J*xdB(x)  , A <_  B , and  <_  S2  , then 


<_  for  all  k = 1,2,  ...  , 


where  D (D  ) is  the  delay  in  queue  of  the  k customer  in  system 

4C  iC 

A/s  /m  (B/S2/m)  . 


Proof: 


Let  T^(T^)  , i = 1,2,  ...,  be  the  interarrival  time  of  the 

i*"*1  and  (i  + l)1"*1  customer  in  system  A/S^/m  (B/S2/m)  . Since 

T1  ± Ti  and  E(T  ) = E(T^)  for  all  i *=  1,2,  ...  , we  have  by 

v 

property  (c)  of  variability  that  -T  <_  -T^  for  all  i * 1,2 

Let  S = s^  , 1*1,2  > with  probability  one.  Noting  t..at  customers 
depart  in  the  order  of  their  arrivals  in  a parallel  server  queue 
with  constant  service  times,  we  have  the  following  relations: 


(3.4.2) 


max  0,D,  + s 

I k-m 


k-i  I 

1-  * tj 

j-k-m 


for  all  k > n , 


Remark: 

The  same  proof  will  go  through  if  we  allow  the  service  times 
to  be  random  but  satisfying  the  following  conditions: 


(i)  S1  £ S2 


(ii)  P(a.  1 si  1 b > « 1 , i * 1,2  , 
(iii)  A(b^  - a^)  = 0 and  B(b2  - 


In  other  words,  we  want  the  customers  to  depart  in  the  order  of 
their  arrivals.  Under  this  condition,  relation  (3.4.2)  and  (3.4.3) 
remain  valid  and  the  inductive  proof  still  holds. 

Finally,  combining  Theorem  3.4.3,  3.4.4,  and  3.4.5,  we  have 
the  following 

Theorem  3.4.6: 

Consider  a tandem  queueing  system  with  n stations  in  series 
where  each  station  can  have  either  one  server  with  general  service 
distribution  or  any  number  of  servers  in  parallel  with  constant 
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service  times.  For  two  such  systems,  if  the  interarrival  and  service 
distributions  are  ordered  in  the  same  way  as  in  Theorem  3.4.3,  3.4.4, 
and  3.4.5,  then  similar  order  relations  hold  for  the  total  waiting 
time  in  system  for  every  customer. 

Proof : 

It  can  be  shown  by  induction  that  the  time  until  any  customer 
leaves  the  whole  system  is  an  increasing  convex  function  of  all  the 
interarrival  and  service  times  of  all  customers  up  to  him.  Hence, 
the  conclusion  follows  in  a similar  fashion  as  in  the  proofs  of 
Theorem  3.4.3,  3.4.4,  and  3.4.5.B 


3,5  Examples 

The  results  obtained  in  the  last  section  are  quite  useful  for 
finding  bounds  for  total  waiting  times  in  tandem  queueing  systems. 
Suppose  we  know  the  ordering  between  two  systems  and  one  of  them 
can  be  analyzed  exactly.  Then,  bounds  for  the  corresponding  quantities 
in  the  other  system  can  be  readily  obtained.  We  give  two  specific 
examples  below. 


I.  Consider  the  following  four  systems  of  tandem  queues: 


(A)  M/A  /I  -*■  A„/l  •+■  •••  A /I 

i n 

(B)  M/B./l  -*■  B-/1  -*••••-*  B /I 

1 £■  n 

(C)  M/C.  /I  + C,/l  C /I 

l n 

(D)  M/D,  /I  -*■  D,/l  -*•  •••  -*■  D /I 

I / n 


where 


A 

B 

C 

D 


i 

i 

i 

i 


, i 
, i 
, i 
, i 


1 j 2 | * • • | n 

1 } 2 | ■ < * | n 

1 1 2 | ■ • 1 1 n 

1 1 2 j ■ • > j n 


have  NWUE  distributions, 
have  exponential  distributions, 
have  NBUE  distributions, 
are  deterministic. 


Assume  that  ^ • D^  have  the  same  mean  for  all 

i - 1,2,  ...,  n and  the  arrival  processes  are  Poisson  with  the  same 

rate  for  all  four  systems.  It  follows  from  property  (e)  of  vari- 

v v v 

ability  ordering  that  Di  ^ A^  for  all  i «=  1,2,  ...,  n . 

Let  Wa^WB,''c,WD^  stationary  total  waiting  time  of  a 

customer  in  system  A(B,C,D)  . By  Theorem  3.4.4,  we  have 


W, 


v 

< W„ 


D - C - 


v 

< W. 


v 

< W 
B - A 


The  point  of  this  example  is  that  the  distributions  of  and 

Wg  are  known  exactly.  Hence,  we  have  both  upper  and  lower  bounds 
for  E^W^  j anc^  a l°wer  hound  for  E for  all  real  number  r >_  1 

II.  Tembe  and  Wolff  showed  that  the  expected  delay  in  front  of 
the  second  server  in  a M/D/1  -+■  M/1  system  is  bounded  above  by 
p'(y  \ x)  (®ee  (2.2.1))  where  X is  the  arrival  rate  and  p is  the 
service  rate  at  the  second  station. 


Theorem  3.4.4  implies  that  the  corresponding  expected  delay, 

* 

d , is  smaller  for  a M/D/1  -*■  G/l  system  where  G has  NBUE  dis- 
tribution and  the  same  arrival  and  service  rates  are  assumed. 
Therefore,  combining  with  Theorem  2.3.5,  we  have 
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* 

d 


2(1  - X/u)  J 


where 


2 

o 

a 


1/X2  and  02 
g 


is  the  variance  of  G . 


3.6  Some  Stronger  Comparisons 

Sometimes  it  is  of  interest  to  know  not  only  comparisons  of  total 

waiting  times  in  system  but  also  the  delays  in  individual  stations 

in  tandem  queues.  Under  stronger  assumptions  (e.g.,  stochastic 

ordering)  between  interarrival  and  service  times  of  two  systems, 

some  results  in  this  direction  can  be  obtained. 

Consider  two  systems  of  tandem  queues  F/G/l  -*■  H/l  and 

F/G/l  ■+•  H/l  where  F(F)  is  the  interarrival  distribution,  H(H)  is 

the  service  distribution  of  the  second  station,  and  G is  the  service 

distribution  of  the  first  station  for  both  systems.  For  system 

F/G/l  -*•  H/l  , let  Dn|Dnj  be  the  delay  of  the  ntl1  customer  in  front  of 

the  first  (second)  station,  S and  R be  the  service  times  of  the  n 

n n 

customer  at  the  first  and  second  station  respectively,  and  be  the 

til  til 

interarrival  time  of  the  n ‘ and  (n  + 1)  customer.  Let  all  barred 
notations  |e.g.,  denote  the  corresponding  quantities  in  system 

F/G/l  -*•  H/l  . 

From  (2.3.8),  we  have  that 


(3.6.1)  D , - max 
n+i 


* 

0,  min  D + S + D 
l n n n 


+ R - T - S . - ,D  + R - S . 
n n n+1  n n n+1 


]]■ 


and 


_* 

(3.6.2)  D = max 
n+1 


[0,  min  D + S + D*  + R - T - 
n n n n n 


_* 


S . . ,D  + R - S . 
n+1  n n n+1 


]] 
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Denote  stochastic  ordering  by  <_  . Before  proving  the  next 
theorem,  we  need  the  following 

Lemma  3.6.1: 

Let  X , Y , i - 1,2,  n , be  independent  random  variables 
st 

Then,  X^  <_  for  all  i - 1,2,  ...,  n if  and  only  if 

f(X1.X2,  V - f(Yl,Y2 V 

for  all  increasing  function  f . 

Proof ; 

Straightforward  by  induction  on  n . ■ 

Theorem  3.6.2; 

st  - st  - 

If  F > F and  H H , then 

* st  _* 

Dfc  <_  for  all  k « 1,2,  ...  . 


Proof : 


Conditioning  on  the  sequence  of  service  times,  i * at 

first  station,  it  follows  from  (3.6.1)  and  (3.6.2)  that  D* ^D*  j 

are  increasing  functions  of  (-T1 -Tk-1’R1’  \ 1^ 

((-T^,  ....  -Tk_^,R^,  ...,  R^_^))  for  all  k - 1,2,  ...  . Hence 
by  Lemma  3.6.1, 

K 1 (Sj’j-l]  [6k  1 { SJ 1 j -l]  ,ot  *n  k ■ !.2 


Unconditioning  the  above  finishes  the  proof.! 
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Theorem  3.6.3: 

If  F 5^  F and  H <_  H , then 


D*  < D*  for  all  k = 1,2, 


Proof : 


Let  f be  an  arbitrary  increasing  convex  function.  Denote 


* Ft1,  ....  f. 


, and  G. 


lT  .....  T * f » and  Gc  c the  joint 

1’  ’ k-1  Al»  1ic_i  si»  •••.  sk 


distribution  of  and  respectively. 

Conditioning  on  and  , „e  have 


e f v 1 ‘vi&vU  • 


*•••  Tk-l(tl’  tk-l>dGS1.  ....  S,(61 sk> 


iffikz) 1 “i’S-vi-i!  • 


(t-  ) • • • ) t.  )dG 


— T T * • • • » r (s.  , • • • • s.  ) 

1*  ik-l  1 k*1  si»  •••*  Sk  1 k 

since  [£(d‘)  | end  [f(B*)  | 

are  increasing  convex  functions  of  (R^ , . . . , and  (R^ , . . . , r^  ) 

respectively  (see  (3.6.1)  and  (3.6.2)) 

-//E!fK)  i (ti>w-<*i>wl  • 


t * • • • $ t,  . )dG  (s  , ....  s ) 

Al»  •••»  1jc_j  1 k-1  S , . . . , S,  1 It 


by  Lemma  3.6.1 


E{f(Dk)}  .1 
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Remark: 

(i)  In  the  proofs  of  Theorem  3.6.2  and  3.6.3,  we  only  needed 

6 1 

(Tlf  •••»  Tk  L (T ^ , . ...  for  all  k ■ 1,2,  ...  . 

Hence  the  arrival  processes  may  be  "dependent."  Also,  only 
st  - — 

(R^ , . . . , R^  (R^ i • « » » for  all  k * 1,2,  ... 

was  used  in  Theorem  3.6.2. 

(ii)  Analogues  of  Theorem  3.6.2  and  3.6.3  can  easily  be  proven 
for  two  systems  with  n (n  >_  2)  queues  in  tandem  where 
the  first  through  (n  - 1)C^  servers  have  the  same  service 
distribution  for  both  systems. 
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CHAPTER  4 

A HETEROGENEOUS  ARRIVAL  QUEUEING  LOSS  MODEL 

4.1  Descriptions  of  the  Model 

Consider  a single  server  queueing  loss  system.  The  system 
alternates  between  two  phases  (1  and  2).  During  phase  i , i - 1,2  , 
customers  arrive  according  to  a Poisson  process  with  rate  X^  , 
and  the  service  times  required  by  customers  have  distribution  G^ 
with  finite  mean.  The  amounts  of  time  the  system  spends  in  phase  i , 
i * 1,2  , are  exponential  random  variables  with  rate  ca^  , J4  0 , 
i • 1,2  . An  arriving  customer  who  finds  the  system  busy  departs 
Immediately  and  is  called  a lost  customer. 

Notice  that  the  average  arrival  rate,  X , of  customers  to  the 
system  is 


X 


Xl°2 


X 

+ a 


which  is  independent  of  c . However,  as  c increases,  the  rate 
at  which  the  system  switches  between  phases  increases,  and  the  arrival 
process  converges  to  a stationary  Poisson  process  with  rate  X as 
c goes  to  infinity. 

Let  L(c)  be  the  long  run  proportion  of  customers  lost  for  a 
given  c in  this  model.  In  the  next  section,  we  will  find  an  explicit 
formula  for  L(c)  , and  examine  the  properties  of  it  in  some  special 


i 


cases . 


A.  2 The  Loss  Formula 


The  system  operates  similar  to  a counter  model.  Each  time 
an  arriving  customer  finds  the  system  empty,  he  enters  service 
immediately  and  block  the  system  for  a random  period  of  time  whose 
distribution  depends  on  which  phase  the  system  is  in  at  the  time  of 
his  arrival.  During  the  blocking  period,  all  arriving  customers 
are  lost.  Hence  the  system  alternates  between  busy  and  idle  periods. 


Let  Z -i,i«l,2,  if  the  system  is  in  phase  i at  the 
beginning  of  the  nC^  busy  period.  We  will  start  by  proving  the  following 


Theorem  A. 2.1: 


The  embedded  Markov  chain  {Z^  , n = 1,2,  ...}  has  a stationary 
distribution  and  is  given  by 


(A. 2.1) 


\ - lta  P!z„  - 11 ' p„  + p„  ■ 
n-*®  12  21 


(A. 2. 2) 


’2  - ^ p<z„  ■ z>  - ?— nr-  . 

n-*°°  12  21 


where 


(A. 2. 3)  P, 


2 

°1*1*2  + C01*2  + cala2*2  ~ °‘l*l*2l*>l^C  ^al  + °2^ 
(c^  + + ca2^1  + 
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(4.2.4)  P 


21 


a2*lX2  + ca2*l  + cctla2Xl  “ a2Xi*2lMc^al  + 
(a1  + c^MA^  + ca2Xl  + ca1X2^ 


and 


t'lUO*!  + a2)]  « j e 


-c(a1+n2)x 


dG1(x)  , 


t>2 1 c Cai  + a2>]  •=  J e 


-c(a1+a2)x 


dG2 (x) 


Proof : 

Let  the  transition  probability  matrix  of  the  Markov  chain 
{Z^  , n = 1,2,  . . .}  be 

?12\ 

\P21  P22  / 

where  P^  = P{Zn+1  - j | Zfl  - i}  , i ■ 1,2  , j - 1,^  , 

If  P^'s  were  given,  then  * (tt^t^)  can  be  obtained  by 
solving  Z anc*  + ^2  = ^ " ii  is  easy  to  verify  that 

(4.2.1)  and  (4.2.2)  indeed  satisfies  the  above.  Hence,  we  only  have 
to  find  P , i ■ 1,2  , j ■ 1,2  . 

Let 


R F * the  probability  that  the  system  will  be  in  phase  £ 

V l 

at  the  end  of  a busy  period,  given  that  it  is  in  phase 
i at  the  beginning  of  that  busy  period,  i * 1,2  , 

i - 1,2  ; 
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and 


the  probability  that  the  system  will  be  in  phase  j 
at  the  beginning  of  next  busy  period,  given  that  it 
is  in  phase  £ at  the  end  of  a busy  period,  £ ■ 1,2  , 


j * 1,2  . 


Conditioning  on  the  phase  of  the  system  at  the  end  of  a busy 
period,  we  have 


2 

(4.2.5)  P - l PR  E • PF  B , i - 1,2  , J - 1,2  . 

ij  £-1  V £ E£’Bj 

We  have  to  find  P and  P , i ' 1,2  , 1 c 1,2  . 

VE£  E£’Bj 

Let  Z(t)  = i , i * 1,2  , if  the  system  is  in  phase  i at  time 
t . Then,  (Z(t),  t ^ 0}  is  a continuous  time  Markov  chain.  Denote 
P^(t)  the  probability  that  the  phase  of  the  system  is  j at  time  t , 
given  that  it  is  in  phase  i at  time  0 . It  may  be  readily  verified 
that 


11 


-c(a  +o  )t  f -c(a  +a_)t"l  f a-  1 

(',  ■ * + L1  - * J • • 

■..<■>  ■ [■  - • [^j  ■ 


P22(t)  = e 


and 


21 


<■>  ■ [*  - •'•'•‘“■i  • yy . 
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satisfies  the  set  of  Kolmogorov  backward  differential  equations 
(see  Ross  [20,  p.  Ill)  for  the  process  (Z(t),t  ^ 0}  . Therefore, 
by  conditioning  on  the  length  of  a busy  period,  we  have 


PB  ,E  “/ Pii(t)dGi(t)  . i - 1.2  , t - 1,2  , 
i * « 


or  explicitly, 


a2  al 

“ 7 + 7 * [c(a,  + a,)]  , 

B^,E^  al  a2  1 1 2 


B1’E2  al  + a2 


77  ' <1  - ‘lIc(al  + °2>J)  • 


^7^7  • II  - *2(C<=1  + a2))l 


1 2 

PR  r = 7 + 7 * ♦9[c(a  + a,)]  . 

^2*^2  ai  + a2  al  + a2  ^ ^ ^ 


Next,  we  will  compute  P_  , £=1,2,  j =1,2.  Since 

VBj 


P_  _ = P{an  arrival  occurs  before  the  system  switches  to  phase  2}  • 1 
E1’B1 


+ P{no  arrival  occurs  before  the  system  switches  to  phase  2}  • P 


E .B 
1’  1 


+ /— -1— W— 2— -V  P 

{1  + C°1  \Xl  + cal  / \X2  + Ca2  / E1’B1 
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Solving  for  P_  „ , we  get 

E1’B1 


P „ X1(X2  + Ca2) 
E1,B1  X1X2  + ca2Xl  + CalX2 

By  symmetry,  we  have 


I, 


Also, 


A«  (A.  ca , ) 

p = - ± i 

£2^2  + ca2Xl  + CalX2 


P = 1 - P 

VB2  rErBl 


CalX2 


A^A2  + ca2A^  + ca^^2 


and 


P * 1 - P 

E2’B1  E2’B2 


Ca2Xl 


X1X2  + Ca2Xl  + CalX2 


Substituting  P and  P , i - 1,2  , j - 1,2  , £ - 1,2  , 

V £ V j 

into  (4.2.5)  gives  (4.2.3)  and  (4.2.4).B 

Let  be  the  number  of  arrivals  in  the  iC^  busy  period, 

i “ 1,2 The  distribution  of  depends  on  the  phase  of  the 

system  at  the  beginning  of  the  i^  busy  period.  The  explicit  form 
of  it  will  be  difficult  to  obtain.  However,  we  have  the  following 
lemmas  which  are  needed  in  the  proof  below. 


r 
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Lemma  4.2.2: 

The  Markov  chain  { (Z  ,K  ) , n - 1,2,  ...}  has  a stationary 

r n 

distribution. 

Proof : 

Let  rik  , i * 1,2  , k * 1,2 be  the  probability  that  the 

' total  number  of  arrivals  in  a busy  period,  which  starts  when  the 

phase  of  the  system  is  in  i , is  k . Then,  it  is  easy  to  check 
that 

P{(Z  ,K  ) = (i ,k) } = nr.  , i = 1,2  , k = 1,2,  ... 
n n l IK 

is  a stationary  distribution  of  { (Z  ,K  ) , n = 1,2,  ...}  . H 

n n 

Lemma  4.2.3; 

If  has  the  distribution 

P { K.^  = k}  e "^L^lk  1I2r2k  * ^ c ^ * ...  , 

then  {K  , n = 1,2,  ...}  is  a stationary  ergodic  process, 
n 

Proof : 

Since  a stationary  distribution  for  the  Markov  chain 

{ (Z  ,K  ) , n = 1,2,  ...}  exists  and  the  chain  is  indecomposable, 
n n 

it  follows  from  Theorem  7.16  of  Breiman  [3]  that  the  process 

{ (Z  ,K  ) , n - 1,2,  ...}  with  the  stationary  initial  distribution 
n n 

given  in  Lemma  4.2.2  is  ergodic.  Therefore,  {K^  , n ■=  1,2,  ...} 

is  also  a stationary  ergodic  process  since  is  just  a coordinate 

mapping  of  (Z  ,K  ) (see  Theorem  6.31  of  Breiman) . ■ 
n n 
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Now  we  are  ready  for  the  main 


Theorem  4.2.4: 


L (c)  - 1 - 


E (K^) 


where 


oo  t 


E(KX)  - 1 + "iff  tXlPH(y)  + x2P12(y)]dydGl(t) 
0 0 


(4.2.6) 


00  t 


*2  j j [^iP2i(y)  + X2P22^y^dydG2^  ' 


0 0 


Proof : 

Observe  that  in  each  busy  period,  only  the  first  customer  can 
enter  the  system  and  receive  service,  the  rest  of  them  are  lost. 
Hence  the  proportion  of  customers  lost  in  n busy  period,  L (c)  » 
is  given  by 


Ln(c) 


l K - n 
1-1  1 
n 

l Ki 

i-l  1 


1 - 


n 

l Ki 

i-l  1 


- 1 - 


l y« 

i-l  1 
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Letting  n goes  to  infinity,  we  get  the  long  run  proportion  of 
customers  lost  which  is  also  the  stationary  probability  that  an 
arriving  customer  will  be  lost 


L(c) 


lim  L (c) 

n-Kc  n 


lim  £ K /n 

n-K»  i«l 


= 1 


1 

Etttp  ' 


The  last  equality  follows  because,  without  loss  of  generality, 

we  can  assume  that  {K  , n - 1,2,  ...}  is  a stationary  ergodic 

n 

process  by  letting  K.^  have  the  stationary  distribution  given  in 
Lemma  A. 2. 3,  and  hence,  by  ergodic  theorem,  (e.g..  Theorem  6.28  of 
[3]) 


lim 

n-H» 


l K< 


i-1 


E(KX) 


W.P.l. 


So,  we  only  have  to  compute  E(K^)  . 

Let  N , i = 1,2  , be  the  total  number  of  customers  who  arrive 
in  a busy  period  with  length  T and  starting  with  the  system  in 
phase  i . For  i ■ 1,2  , j * 1,2  , define 


(y) 


if  the  phase  of  the  system  at  time 
j , given  that  at  time  0 it  is  i 

otherwise, 


y 


is 
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then  T m f I..(y)dy  is  the  amount  of  time  the  system  spends 
ij 

in  phase  j in  [0,T]  if  it  starts  from  phase  i . Therefore, 


E(Nt) 


XlE(Til)  + X2E(Ti2) 


'T 

■ V / 

.0 


In(y)dy  + X2E  I Ii2(y)dy 


X X 

X-j J P A1 (y)dy  + x2%/*Pi2(y)dy  * 


Hence  (4.2.6)  follows  by  conditioning  on  the  length  and  the 
phase  of  the  system  at  the  beginning  of  a busy  period.® 


4.3  Special  Cases 


Case  1:  c * 


Let  G,  and  G_  have  the  same  mean  — . As  we  mentioned 

12  p 

earlier  in  Section  4.1,  the  parameter  c in  this  model  controls 
the  rate  at  which  the  system  switches  between  phases.  And,  as  c 
goes  to  infinity,  the  arrival  process  converge  to  a stationary 
Poisson  process  with  rate  X . Indeed,  from  (4.2,6),  we  have 


lim  E(K, ) -*•  1 + - 
1 u 


or  equivalently, 


lim  L(c) 


X + u 
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which  is  exactly  the  loss  formula  for  a no  queue  allowed  M/G/l 
loss  system  with  arrival  rate  X and  service  rate  y . 


Case  2:  X2  * 0 

The  expression  we  obtained  for  L(c)  is,  unfortunately,  too 
complicated  to  determine  whether  it  is  a decreasing  function  of  c . 
Some  numerical  experiences  indicate  that  L(c)  is  decreasing  for 
Erlang,  exponential  and  hyperexponential  service  distributions. 


One  interesting  special  case  is  to  let  X,,  ■=  0 . Note  that, 
in  this  case,  the  arrival  process  becomes  a renewal  process,  and  is 
sometimes  called  an  interrupted  Poisson  process  in  teletraffic 
terminology.  It  arises  naturally  as  an  over  flow  process  from  some 
queueing  loss  systems.  Letting  X2  - 0 in  (4.2.6),  we  have 


(4.3.1)  E(K:) 


1 + 


/ Vi 

lc(ai  + a2) 


[1  - ♦ [c(a.  + a,)]]  + ± 

112  y 


which  may  be  shown  by  differentiating  to  be  a decreasing  function 
of  c . Thus,  L(c)  is  also  a decreasing  function  of  c . Hence, 
in  this  system,  with  a given  service  distribution  , the  more 
stationary  the  arrival  process  is  the  smaller  L(c)  will  be. 

For  a fixed  c , we  can  also  compare  L(c)  for  two  such 
systems  with  different  service  distributions  and  H2  . Suppose 

and  H2  have  the  same  mean  and 

00  00 

<4.3.2)  J e‘stdH1(t)  < Je"8tdH2(t)  , 

0 0 
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i.e.,  the  Laplace  transforms  of  and  are  ordered.  Then, 

(4.3.1)  and  (4.3.2)  imply  that 

L1(c)  > L2(c)  , 

where  L^(c)  , i-1,2  , is  the  loss  formula  for  the  system  with 
service  distribution  . 

v 

One  sufficient  condition  for  (4.3.2)  to  hold  is  <_ 
and  J*xdH^  (x)  since  e St  is  a convex  function  of  t . 

Therefore,  for  interrupted  Poisson  arrival  process,  the  more  regular 
the  service  distribution  is  the  "worse"  the  system  performance  will 
be  for  a single  server  loss  system.  This  is  perhaps  a surprising 
phenomenon  and  we  will  look  into  this  further  in  the  next  section. 

4.4  Some  Related  Results 

Regularity  of  arrival  process  and  service  times  seems  to  work 
to  good  effect  with  respect  to  the  usual  performance  measures  in 
a number  of  queueing  systems.  Therefore,  it  is  natural  to  ask: 

Does  this  general  principle  holds  for  loss  systems? 

We  consider  two  types  of  comparisons.  One  can  either  compare 
two  systems  with  different  service  distributions  for  fixed  arrival 
process,  or,  for  fixed  service  distribution,  compare  two  systems 
with  different  interarrival  distributions.  Counterexamples  of  both 
types  will  be  presented  below.  Hence,  the  answer  to  the  above 
question  is:  no. 

Consider  a single  server  loss  system.  All  arriving  customers 
who  find  the  system  empty  are  lost.  Let  the  service  times  be 
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deterministic,  say  10  minutes  for  each  customer.  Compare  two  such 
systems:  one  has  i.i.d.  interarrival  times  T^  , i * 1.2,  ...  , 
such  that 

7 minutes  w.p.  1/2 
11  minutes  w.p.  1/2  , 

and  one  with  deterministic  interarrival  times  with  = 9 minutes 
for  all  i = 1,2,  ...  . Then,  it  is  easy  to  see  that 

Proportion  of  customers  lost  in  first  system  « 1/3 
Proportion  of  customers  lost  in  second  system  *=  1/2  . 

Thus,  the  second  system  which  has  more  regular  interarrival  times 
behaves  "worse”  than  the  first  system. 

Another  example  in  the  same  vein  is  obtained  by  reversing 
the  roles  of  interarrival  and  service  time  distributions  in  the  above. 
Consider  two  loss  systems  with  deterministic  interarrival  times, 
say  10  minutes  for  each  customer  to  arrive.  Let  the  service  times 
, i = 1,2,  ...  for  customers  in  the  first  system  be  such  that 

( 9 minutes  w.p.  1/2 

Si  = 

(13  minutes  w.p.  1/2  , 

and  = 11  for  the  second  system.  Again,  we  have 

Proportion  of  customers  lost  in  first  system  ■ 1/3 
Proportion  of  customers  lost  in  second  system  * 1/2  . 


J 


Hence,  the  system  with  more  regular  service  times  has  a "larger" 
loss.  Note  that  the  system  discussed  at  the  end  of  the  last  section 
also  exhibits  the  same  contrary  behavior,  but  the  example  here 
is  simpler. 

The  next  question  is:  Is  there  any  condition  under  which  a 
loss  system  with  more  regular  interarrival  and  service  distributions 
does  perform  better?  Some  partial  answers  to  this  question  for 
single  server  loss  systems  are  given  below. 

Consider  two  single  server  loss  systems  F^/G/l  and  F2/G/I 
where  F^  , F2  are  interarrival  distributions  and  G is  the  service 
distribution  for  both  systems.  Let  L^(L2)  be  the  proportion  of 
customers  lost  in  system  F^/G/l  (Fj/G/l)  . 

Theorem  4.4.1: 


Suppose  / xdF  (x)  =/xdF_(x)  = 1/A  , F (F,)  has  NBUE  (NWUE) 

0^0  1 
00 

distribution,  and  f xdG(x)  K 1/p  . Then, 

0 


L1  - A + p - L2 


where  - +-  ^ is  the  corresponding  loss  formula  for  a M/G/l  loss 
system  with  arrival  rate  A and  service  rate  p . 


Proof: 


Let  m^Ct)  be  the  expected  number  of  arrivals  in  [0,t]  for 


the  i system,  i * 1,2  . From  Theorem  4.2.4,  we  have 
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Li  ■ 1 - EtcTy 


. i - 1,2 


where  E(C  ) » 1 + f m (t)dG(t)  is  the  expected  number  of  arrivals 
0 

in  a busy  period. 

Hence,  we  only  have  to  show 


(4.4.1)  J * m^  (t)dG(t)  <_  Xt  <_J  m2  (t)dG(t)  . 

0 0 

But,  (4.4,1)  follows  directly  from  a well-known  fact  from  renewal 
theory  that  m^t)  _<  Xt  <^m2(t)  . ■ 

Theorem  4.4.2: 


Suppose  /e  dF  (x)  <_/e~  dF,(x)  , / xdF.  (x)  =/xdF0(x)  , 
0 -L0  0 1 0 

and  G has  hyperexponential  distribution.  Then, 


^1  — l2 


Proof : 


Again,  we  only  have  to  prove 


w w 

J m1(t)dG(t)  £ j m2 (t)dG(t)  . 


Let  $ (x)  * J e sxdF.  (x)  , i - 1,2  . Now,  since  m (t)  - \ F[n^(t)  , 

0 1 n-1  1 


we  have 
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“i(s)  = • 1 * X'2  ’ 


where  m^(s)  is  the  Laplace  transform  of  m^Ct)  . Therefore, 
$^(s)  <_  4>2(s)  ++  n^(s)  £m2(s)  . Integrating  by  part,  we  get 


m1(s)  e Stdm(t) 


/"  i 


(t)se  Stdt 


Thus,  f m (t)se  Stdt  <_ f m (t)6e_stdt  for  all 
0 1 0 i 


s . 


-St 


Notice  that  se  is  exactly  the  density  of  an  exponential 
distribution  with  rate  s . It  follows  that 


j m1(t)dG(t)  <_J  m2(t)dG(t) 


if  G is  a mixture  of  exponential  distributions. 
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